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Abstract 

The production of Kg mesons and A baryons in quark and gluon jets has been 
investigated using two complementary techniques. In the first approach, which 
provides high statistical accuracy, jets were selected using different jet finding al- 
gorithms and ordered according to their energy. Production rates were determined 
taking into account the dependences of quark and gluon compositions as a function 
of jet energy as predicted by Monte Carlo models. Selecting three-jet events with 
the k± (Durham) jet finder (j/cut = 0.005), the ratios of Kg and A production rates 
in gluon and quark jets relative to the mean charged particle multiplicity were found 
to be 1.10 ± 0.02 ± 0.02 and 1.41 ± 0.04 ± 0.04, respectively, where the first uncer- 
tainty is statistical and the second is systematic. In the second approach, a new 
method of identifying quark jets based on the collimation of energy flow around the 
jet axis is introduced and was used to anti-tag gluon jets in symmetric (Y-shaped) 
three-jet events. Using the cone jet finding algorithm with a cone size of 30°, the 
ratios of relative production rates in gluon and quark jets were determined to be 
0.94 ± 0.07 ± 0.07 for Kg and 1.18 ± 0.10 ± 0.17 for A. The results of both analyses 
are compared to the predictions of Monte Carlo models. 
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1 Introduction 



In recent years, clear differences between quark and gluon jets liave been establislied 
experimentally These studies have mostly exploited the large samples of TP decay 

events recorded at the CERN LEP collider, and have examined three-jet events with a 
symmetric event topology favouring the direct comparison of quark and gluon jets of the 
same energy produced in the same event topology. In particular, gluon jets have been 
measured to have a larger mean particle multiplicity, a softer fragmentation function 
and a larger angular width than quark jets of the same energy. Comparisons between 
jets produced in e+e~ and pp collisions also indicate that gluon jets are broader than 
quark jets [0. This is in qualitative agreement with the predictions of perturbative QCD 
1^. Perturbative QCD makes no explicit predictions for the individual hadron species 
produced in these jets, but some QCD Monte Carlo models of hadronization predict that 
the relative hadron production rates in quark and gluon jets differ for different hadron 
species. 

Few experimental results are available on the production of identified particles in quark 
and gluon jets. Results from e"'"e~ annihilations at center-of-mass (cm.) energies around 
the T(15') resonance (~ 10 GeV) indicate that baryons are produced about 2.5 times more 
copiously in direct T(IS') decays (T(15') — >• ggg hadrons) than in continuum events 
(e"'"e~ — >■ qq — >■ hadrons) whilst no such enhancement was observed for mesons [^j. An 
OPAL investigation |]T0| studied identified particle production in jet topologies enriched 
in gluon and quark jets and found no jet-dependent differences in the production of 
mesons and charged particles (other than protons). In contrast, baryons were found to be 
produced more copiously in gluon-enriched jet samples. The DELPHI collaboration has 
reported measurements of K^, K°, p and A in secondary- vertex tagged quark and gluon 



jets in symmetric (Y-shaped) events [|r^ . Within relatively large statistical and systematic 
uncertainties, they find the ratio of identified particle rates in quark and gluon jets to be 
consistent with that for charged particles. Recently, L3 has concluded that production of 
K° and A in both quark and gluon jets is well modelled by string fragmentation [ll2| . 

In this paper an experimental comparison of K° and A production in quark and gluon 
jets is presented. Two complementary analyses of data recorded with the OPAL detector 
at LEP are presented, using different approaches to identify quark and gluon jets. One 
analysis separates quark and gluon jets according to their energies, giving a large sample 
of events, whilst the other selects a smaller sample of tagged quark and gluon jets in 
symmetric (Y-shaped) events, which allows for a simpler interpretation, but with larger 
uncertainties. 

In the first analysis (the 'energy-based analysis') the production of Kg mesons and A 
baryonsQin quark and gluon jets was investigated in three-jet events of different topologies 
selected with the fc_L (Durham) |13| or a cone jet finder. The jets were ordered by their 



energy since the lowest energy jets are mainly induced by gluons, and higher energy jets 
mainly by quarks. Motivated by Monte Carlo investigations a similar energy dependence 
for the production rates of all hadron species was assumed, and the production rates of Kg 



^For simplicity A refers to both A and A. 
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and A relative to those of charged particles were determined. The experimental relative 
rates were corrected for the underlying mixture of quark and gluon jets, allowing the Kg 
and A relative production rates in pure quark and gluon jets to be obtained. 

In the second analysis (the 'Y-event analysis'), a comparison was made of the absolute 
production rates of Kg mesons and A baryons in quark and gluon jets produced under 
the same conditions, embedded in similar event topologies. Symmetric three-jet events 
were analysed, where the two lower energy jets (assumed to be initiated by a quark and 
a gluon) were produced at about 150° with respect to the higher energy jet. A sample of 
anti-tagged gluon jets containing about 30% of the symmetric event sample was isolated 
by means of a new method of identifying the quark jets based on the observation that light 
quark jets are more coUimated than gluon jets 0. The inclusive yield of charged particles 
in these jets was also measured, allowing relative rates to be evaluated. The production 
rates of the particles in the lower energy jets of the inclusive symmetric sample were 
determined, and a correction applied in order to obtain measurements corresponding to 
pure samples of quark and gluon jets. 

Section ^ gives a description of the main features of the OPAL detector, of the data 
and simulated event samples and of the reconstruction algorithms for Kg and A. Section |] 
describes the details of the energy-based analysis, and section H the Y-event analysis. The 
results of both analyses are presented and discussed in section |^, before the summary is 
given in section 



2 The OPAL Detector and Data Samples 
2.1 The OPAL Detector 



The OPAL detector is described in detail elsewhere Of most relevance for the present 
analyses are the tracking system and the electromagnetic calorimeter. The tracking sys- 
tem consists of a silicon microvertex detector, an inner vertex chamber, a large-volume 
jet chamber and specialised chambers at the outer radius of the jet chamber which im- 
prove the measurements in the z direction| (z-chambers). The tracking system covers 
the region | cos^| < 0.95 and is enclosed by a solenoidal magnet coil with an axial field 
of 0.435 T. The tracking detectors provide momentum measurements of charged particles, 
and particle identification from measurements of the ionisation energy loss, dEjdx. Elec- 
tromagnetic energy is measured by a lead-glass calorimeter located outside the magnet 
coil, which covers I cos^l < 0.98. 



2.2 Data Samples 

The energy-based analysis, which is not statistics limited, used 2.8 million events recorded 
between 1992 and 1994 whilst the Y-event analysis used the full OPAL data sample of 

^The coordinate system is defined so that z is the coordinate parallel to the e~ beam axis, r is the 
coordinate normal to the beam axis, </> is the azimuthal angle around the beam axis and Q is the polar 
angle with respect to z. 
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about 4.2 million hadronic events collected around the Z" peak from 1990 to 1995. The 
procedures for identifying hadronic events using measurements of tracks and electromag- 
netic energy are described in |TB[. The criteria applied to select tracks and deposits of 
electromagnetic energy (clusters) for the analyses were identical to those in 0. Each 
accepted track and unassociated electromagnetic cluster was considered to be a particle. 
Tracks were assigned the pion mass and electromagnetic clusters were assigned zero mass 
since they originate mostly from photons. 

To reduce background from non-hadronic decays of the Z° and to eliminate events in 
which a significant number of particles were lost near the beam direction, both analyses 
required that | cos(6'thrust)| < 0.9 (6'thrust is the polar angle of the thrust axis); also, there 
had to be at least five accepted tracks. The residual background from all sources was 
estimated to be less than 1%. In addition, for the energy-based analysis, events were 
rejected if they contained tracks with measured momentum greater than 60 GeV/c, if the 
absolute value of the vector sum of all selected particles |ptot| exceeded 30 GeV/c, or if 
the visible energy (the sum of the energies of all the accepted tracks and clusters) was 
less than 40% of the center-of-mass energy. 



2.3 Reconstruction of Kg and A 

The neutral strange Kg mesons and A baryons were reconstructed by their decay channels 
Kg — s> 7r+7r~ and A — > vr^p. The reconstruction algorithms, signal definitions, efficiency 
corrections and background subtractions are all described in [16| and [1^. Briefly, tracks 
of opposite charge were paired and regarded as a secondary vertex candidate if at least one 
track pair intersection in the plane perpendicular to the beam axis satisfied the criteria 
of a neutral two-body decay with the appropriate lifetime. Each track pair passing these 
requirements was refitted with the constraint that the tracks originated from a common 
vertex, and background from photon conversions was suppressed. For A candidates, in- 
formation from dE/dx measurements was used to help identify the vr and p for further 
background suppression. Two sets of cuts are described in |jl^ for A identification.^ For 
the energy-based analysis, A candidates were reconstructed using method 1, which is op- 
timized to have good mass and momentum resolution, while in the Y-event analysis the 
more efficient method 2 was employed to maximize the number of A candidates. 

Candidates for Kg with momentum greater than 0.150 GeV/c and with an invariant 
mass in the range 0.3 GeV/c^ < 772,7^,7 < 0.8 GeV/c^ and A candidates with momentum 
greater than 0.520 GeV/c and with an invariant mass in the range 1.08 GeV/c^ mT^p < 
1.2 GeV/c^ were retained for further analysis. The Kg rates were determined by fitting the 
mass spectrum with a third-order polynomial excluding a signal region of ±0.05 GeV/c^ 
around the nominal mass. The shape of the background under the A signal was fitted 
using a function of the form (1 — e-«(mTp-i.077)-j ^ . ^^^^ _j_ excluding a signal window 
of ±0.012 GeV/c^ around the nominal A mass. For each particle species the entries in the 
signal region were summed and the background as determined by the fit was subtracted. 
The efficiency of the identification algorithms was determined as a function of candidate 



^The two sets of cuts are described as 'raethod 1' and 'method 2' in 111 
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momentum from Monte Carlo events, and subsequently used to correct the number of 
observed signal events. 

For the energy-based analysis, candidates for A and Kg were formed and then supplied 
to the jet finder, whereas this order was reversed in the case of the Y-events. Monte Carlo 
studies showed that for the Y-event analysis there is no systematic effect due to the order 
of jet finding and particle reconstruction. 



2.4 Monte Carlo Event Samples 



Samples of Monte Carlo hadronic events with a full simulation of the OPAL detector [|I| 
and including initial state photon radiation were used to evaluate the detector acceptance 
and resolution, and to study the efficiency and purity of the quark and gluon jet identi- 
fication and the particle reconstruction algorithms. In total, 7 million simulated events 



were available, of which 4 million were generated by Jetset 7.4 |[T^ with fragmentation 
parameters described in [Q, and 3 million generated by Jetset 7.3 with fragmentation 
parameters described in 0. The Jetset 7.4 events included updated particle decay ta- 
bles and heavy meson resonances and were processed using a more recent version of the 
detector simulation compared to the Jetset 7.3 sample. It should be noted that there 
are significant differences in the simulation of baryon production between the two Jetset 
samples. 

For comparison with the experimental results, the Monte Carlo models Jetset 7.4 
and Herwig 5.90 were used. The models both give a good description of global event 



shapes and many inclusive particle production rates, but differ in their description of the 
perturbative phase and their implementation of the hadronization mechanism. 

Tracks and clusters were selected in the Monte Carlo events, which include detector 
simulation, in the same way as for the data: the resulting four-vectors of particles are 
referred to as being at the 'detector level'. Alternatively, Monte Carlo samples with- 
out initial-state photon radiation or detector simulation were used, with all charged and 
neutral particles with mean lifetimes greater than 3 • 10^^*^ s treated as stable. The 
four-vectors of the resulting particles are referred to as being at 'generator level'. The 
remaining quarks and gluons after the termination of the parton shower in these events 
are referred to as being at the 'parton level'. 



3 The Energy-Based Analysis 
3.1 Selection of Three- Jet Events 

Two different types of algorithms are commonly used for jet definition: recombination jet 
finders and cone jet finders. The treatment of high momentum particles which lie close to 

''The fragmentation parameters of Herwig 5.9 were identical to those used by the OPAL tuned version 
of Herwig 5.8 (2l| with the exception of the parameter CLMAX which was set to 3.75 in order to improve 
the description of the mean charged particle multiplicity value in inclusive hadronic Z" decays. 
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the jet axes is similar for both types of jet finder, but there are substantial differences in 
the assignment to the jets of soft particles far from the jet axes. Recombination algorithms 
combine all particles into jets, whereas cone finders do not associate particles outside the 
cones. The recombination jet finders iteratively combine pairs of particles until the scaled 
resolution parameter (y) of all subsequent pairings exceeds the jet resolution parameter 
I/cut- The cone jet finder associates particles into jets that lie within a cone of fixed half- 
angle R. The cone axis is determined from the vector sum of the momenta of the particles 
contained therein. 

To illustrate the sensitivity of the analysis to event topology and jet reconstruction 
algorithm, we selected three samples of three-jet events with different average topologies, 
using either the k± jet finding algorithm [jl3|, or the cone algorithm 0. The main selection, 
containing the largest number of three-jet events, was obtained using the k± algorithm 
with a resolution parameter, ?/cut = 0.005. The second sample of three-jet events (t/win 
sample) was selected from a window of ^/-values between 0.008 and 0.016, i.e., the event 
should be classified as three-jet for some value of ycut within this window. Monte Carlo 
studies at the generator and parton level showed that hadronization effects on the angles 
of jets selected in this sample are small. The average resolution parameter ?/cut for these 
events to be clustered into three jets was about 0.01. Finally, a third sample of three-jet 
events was selected using the cone jet finder with the parameters R = 0.7 rad and e = 
7 GeV, where e is the minimum energy contained in the jet cone. These parameters were 
chosen to give a good correspondence between jets reconstructed at the generator and 
parton levels in Monte Carlo events 0. 

Both jet finding algorithms were applied using reconstructed Kg and A candidates, 
accepted charged tracks (excluding the decay products of Kg and A) and the electromag- 
netic clusters not associated to tracks. Additional quality cuts were then applied to the 
selected three-jet events. Each jet was required to contain at least two charged tracks, 
to have more than 5 GeV of visible energy and to lie in the polar region | cos6'| < 0.9. 
The sum of the angles between all three jets had to exceed 358° to eliminate non-planar 
events, and the angle between the two lower energy jets was required to be larger than 
30°. Finally, the jets in each event were each assigned a calculated energy based on the 
measured jet directions and assuming massless kinematics. 

After all cuts, 24.0% of events were selected by the ?/cut selection, 18.1% by the y^^^ 
selection, and 15.5% by the cone selection. The tighter selection criteria for the ?/„in and 
the cone samples, yielding a lower number of three-jet events, were partly compensated 
by an improved reconstruction quality e.g., a better agreement of the jets at the parton 
and generator level. In total, samples of about 500,000 three-jet events were retained for 
further analysis. In figure |l] the energy spectra of the jets in the ycnt sample is shown. 
The jet energy distributions are more peaked in the ?/cut selection compared to the ?/win 
and cone selections. The relatively small number of events containing a jet with energy 
below 7.5 GeV (~ 15%) were excluded from further analysis. 

In general, the angular separation of the jets (jet topology), is dependent on the event 
selection. In particular, the ?/cut selection gave the most collimated sample of events, with 
an average angle between the two lowest energy jets of about 62°. This can be compared 



9 



to average angles of about 70° for the j/wm sample and about 76° for the cone selection. 



3.2 Determination of Jet Purities 

Since the jets selected in any given energy interval were a mixture of quark and gluon jets, 
Jetset Monte Carlo events were used to determine the quark and gluon jet content of 
the samples in the following manner. The Monte Carlo events were selected in the same 
way as the data events and were accepted if classified as three-jet events at detector level. 
From the partons of these events exactly three jets were reconstructed. The detector level 
jets closest in angle to the parton level jets containing the primary quark and anti-quark 
were considered to be the quark jets and the remaining jet the gluon jet. The term quark 
(gluon) 'jet purity' is defined to be the fraction of jets in a given energy bin, initiated by 
quarks (gluons). The purities of the jets as a function of jet energy are shown in figure |^. 
The lowest energy jet samples contained in excess of 85% gluons, whereas the high energy 
samples were composed of over 95% quark jets. The purity distributions as a function 
of jet energy determined with the k±_ algorithm in the ?/cut and ?/win samples agree to 
within 1%. The gluon jet purity in the cone sample was lower than that of the ?/cut and 
?/win samples by a few percent in very low and in high energy jets. The Herwig Monte 



Carlo model and Jetset in matrix element mode |]19| predict similar jet purities to the 
standard Jetset samples, and detector effects on the purity distribution were found to 
be small. 



3.3 Determination of Particle Rates 

With the Ucut selection, about 100,000 Kg and about 30,000 A were reconstructed. The 
production rates per jet of charged particles. Kg mesons and A baryons in jets, rich, 
rij^o, and ua, are shown as functions of the jet energy in figure |^ after corrections for 
reconstruction efficiency (about 16%) and detector acceptance. The decay products of 
Kg and A were not counted as charged particles. The predictions of the Jetset 7.4 and 
Herwig 5.9 Monte Carlo models are also shown in figure |^. Whereas the Jetset 7.4 
generator describes the experimental data fairly well, the predictions of Herwig 5.9 agree 
poorly with the data. Similar results were obtained with the three-jet events from the 
?/win and the cone selections. 

3.4 Relative Rates and Correction for Jet Impurity 

The particle production rates rise with the jet energy partly due to the changing mixture 
of gluon and quark jets, and partly due to the increased energy available for particle 
production. In order to measure differences due to quark and gluon jets, it is necessary 
to remove this jet energy dependence. If a similar energy dependence for the production 
rates of all hadron species is assumed, then the relative particle production rate (defined 
as the rate of Kg or A production divided by the rate of charged particle production, 
Ry<) = riy^o/rich and Ra = n\/nch) would be independent of the jet energy. Studies 
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of Jetset events showed that there is indeed only a weak energy dependence of these 
relative production rates in pure quark and gluon jets, figure The lines are fits of 
straight lines and have slopes smaller than 2 x 10~^ GeV~^. For gluon jets the slopes 
are smaller than for quark jets, and for Kg smaller than for A. As Jetset gives a good 
overall description of the data over a large range of cm. and jet energies, it was assumed 
for further analysis that the relative particle production rates are constant in the data 
and we choose this assumption to reduce the dependence on Monte Carlo models. The 
consistency of the results with this assumption will be shown later. Particle production in 
jets depends not only on the jet energy but also on the proximity of the other jets in the 
event [^. Further Monte Carlo studies showed that this angular dependence is almost 
the same for all hadron species and does not affect the relative rates. 

The relative particle production rates of Kg and A were computed in each jet energy 
interval from the ratios of the rates shown in figure |^. The relative rates for pure quark and 
gluon jets were unfolded using the jet purities for each three-jet event sample, obtained 
from Monte Carlo results shown in figure ^ A fit was performed to the observed relative 
rates as a function of jet energy Ej using the function: 

RtiE,) = i?J ■ p,iE,) + i?J ■ (1 - p,(E,)), (1) 

where R^{Ej) is the measured relative particle rate, Pg{Ej) the gluon purity of the jets, 
and Rq and Rg the relative rates in pure quark and gluon jets, with h = Kg, A [h = 
hadron). R^ and Rg were assumed to be constant for the reasons stated above. 

The fit was performed in the jet energy range from 7.5 GeV to 45.0 GeV. The relative 
particle rates as a function of the jet energy are shown in figure |^ for the |/cut selection, 
with the lines giving the fit results. The fit results and x^/d-o.f. values for all three 
selections are given in table In order to compare particle production in quark and 
gluon jets, the ratios {Rg/Rq)A and {Rg/Rq)^ were studied. These are given in table |] 
and will be discussed fully in section ^. 



3.5 Systematic Uncertainties 

The following sources of systematic uncertainty on the measured ratios of relative rates 
have been studied. For each source of uncertainty, the difference with respect to the stan- 
dard analysis was used to estimate a symmetric systematic uncertainty. The uncertainties 
listed in table |^ were added in quadrature to arrive at a total systematic error. 

Three-jet event reconstruction: The analysis was repeated with the following changes: 

• charged tracks only were used for the reconstruction of jets, instead of charged 
tracks and unassociated electromagnetic clusters. This check disregards all 
calorimeter information when determining energy flows in the events, so the 
changes represent an extreme situation. A systematic error was therefore de- 
termined from the difference divided by \/V2/2 (= 1.7); 
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the resolution parameter of the k± jet-finder was varied by replacing the ?/cut 
= 0.005 selection by selecting three-jet events in a window of y-values (0.004, 
0.008). The T/win selection was modified by requiring a cut at a fixed value, 

Vent =0.01. 



A and Kg reconstruction: As in [|T6[ and [jl^ the major sources of systematic uncer- 
tainties for the reconstruction of Kg and A were found to be the background de- 
termination and the reliability of the Monte Carlo simulation for distributions on 
which cuts were placed. 



The Kg and A selection criteria were varied as in and |T^. In particular 



the cut on the distance between the reconstructed secondary vertex and the 
first measured hit of the decay particles was loosened from 3 to 9 cm; 

a sideband method [jl^ was applied to determine the backgrounds under the 
Kg and A signals; and 

systematic uncertainties on the strange particle reconstruction efficiencies were 
estimated by calculating the efficiencies separately using the Jetset 7.3 and 
Jetset 7.4 samples (the standard analysis used the combined Jetset 7.3 and 
7.4 samples). The simulation of baryon production differs considerably in these 
samples and they were taken to represent alternative possibilities of generator 
tuning and detector simulation. Studies of events generated using Herwig as 
input to the detector simulation gave efficiencies lying in the range spanned by 
Jetset 7.3 and Jetset 7.4. A symmetric systematic uncertainty was assigned 
using the full difference between Jetset 7.3 and Jetset 7.4 divided by \/l2. 



Quark and gluon jet unfolding: The systematic uncertainties in the determination 
of Kg and A production in quark and gluon jets were obtained by the following 
variations: 



• the upper and lower bounds of the fit range (7.5 — 45.0 GeV) were changed to 
12.5 and 40 GeV respectively; 

• the infiuence of the jet purity determination was studied by varying the purities 
by their systematic uncertainties (about 5%), which were derived as described 
in Briefiy, the uncertainty in the identification of quark and gluon jets was 
estimated by comparing different fragmentation models (Jetset and Herwig) 
and studying detector resolution effects; 

• as a check of the fit procedure, the relative production rates in pure gluon jets 
were obtained by a linear fit to the relative rates rimiPg) as a function of the 
gluon jet purity pg above pg > 8%, and extrapolating to 100% gluon jet purity; 

• the fit function was modified to account for a possible linear energy dependence 
of the relative particle production rates in pure quark and gluon jets, using 
the Jetset slopes as shown in figure It can be seen from table ^ that 
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considering non-zero values of the slopes results in only a small contribution to 
the systematic uncertainties. This means that the relative rates found in the 
experimental data are consistent with the assumption of independence of the 
jet energy. 

The main contributions to the systematic uncertainties came from the variation of the 
cuts on Kg and A candidates, from the differences in detection efficiencies determined 
using Jetset 7.3 and Jetset 7.4, and from the variation of the fit ranges. The ratios 
of relative rates determined from the fit were essentially unchanged if contributions from 
Kg and A decay products were included in the charged particle rates. 

Finally, Monte Carlo events with full detector simulation were analysed in the same 
manner as the data. The relative rates derived from the fits (0.95±0.02 for Kg and 
1.27±0.04 for A) were compared with those determined directly at the generator level in 
these events (0.94 and 1.26 respectively). The agreement for both A and Kg from the ?/cut 
sample is good. The results from the other two jet samples agreed equally well. 



4 The Y-event Analysis 
4.1 Three-jet Event Selection 

For the Y-event analysis, jets were defined with the cone jet finding algorithm, supplying 
all accepted particles as input. The resolution parameters chosen for the jet finder were a 
cone size R = 30°, as in an earlier publication [Q, and a minimum jet energy e computed 
once for each event according to e = 5 x i^vis/ ^/s GeV where ^/s is the cm. energy and 
ii^vis the sum of the energy of the particles. 

The criteria given in references |]l]-§] were followed to select a sample of symmetric 
three-jet events. Each jet was required to contain at least two particles and to lie in the 
polar angle region | cos^| < 0.9, and the sum of the angles between the three jets was 
required to exceed 358°. As in the energy-based approach, the jets in each event were 
assigned a calculated energy based on the interjet angles, assuming massless kinematics. 
Strongly symmetric three-jet events were selected by projecting the jets into the three-jet 
event plane and requiring the angles between the jet with the highest calculated energy 
and each of the two others to be in the range 150° ± 10°. The event plane is defined 
as the plane perpendicular to the sphericity eigenvector associated with the smallest 
eigenvalue. In total, 70,738 symmetric three-jet events were found. The mean calculated 
jet energies were 42.50 ± 0.01 GeV for the highest energy jet and 24.37 ± 0.01 GeV for 
the two lower energy jets. The highest energy jets are likely to be quark jets with high 
probability, due to the nature of the gluon radiation spectrum. From a Monte Carlo study, 
this probability was estimated to be 97.3%. The two lower energy jets were therefore 
assumed to be a quark jet and a gluon jet of equal energy with the same angles with 
respect to the other two jets in the event. The inclusive sample of lower energy jets is 
referred to as the 'normal-mixture' sample of jets. 
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4.2 Gluon Jet Identification 



A gluon-jet enriched sample of the lower energy jets was selected by a new variant of 
the method [0-^ of identifying quark jets in order to 'anti-tag' the gluon jets. For each 
lower energy jet identified as a quark jet, the other lower energy jet was assumed to be a 
gluon jet. This anti-tagged sample of gluon jets was therefore essentially unbiased by the 
tagging method. It contained a well known fraction of gluon jets established from studies 
of simulated events. 

Previous studies employed the reconstruction of secondary vertices or the identification 
of energetic leptons to identify jets that originated from heavy quarks. Whilst this yielded 
an anti-tagged jet sample with a gluon purity typically in excess of 90%, the efficiency to 
identify quark jets was only of order 5%. The present analysis introduces a new method 
to isolate a large sample of quark jets by identifying jets that are collimated, i.e. that 
have a large fraction of their energy close to the jet axis. Studies of quark and gluon jets 
in symmetric events have shown that gluon jets are broader than quark jets 0] and this 
is well reproduced by Monte Carlo models. 

In particular, the quantity f@ is determined, which is given by the fraction of the 
jet's energy contained in a cone co-axial with the jet axis and with half-angle G. If the 
energy of the particles contained in the sub-cone is Eq and the visible jet energy is Ej^t 
then /e = The /e distributions were studied in detector level Monte Carlo events. 
Three-jet events were selected as above, and /e determined for the two lower energy jets. 
Each simulated hadron jet was associated with an underlying quark or gluon jet using 
the method described in Briefly, the two hadron jets that were closest in angle to 
the directions of the primary quark and anti-quark which had evolved from the decay 
were considered to be the quark jets, and the remaining jet was identified as the gluon 
jet. Figure ^ shows the distributions of /e for 6 = 7° for jets in the simulated events that 
were classified as quark or gluon jets. The quark jets tend to have larger /e values than 
the gluon jets, and a sample of jets with a high quark content can therefore be selected 
by requiring /e to exceed some threshold /g,"*, - for example /q"* = 0.75. 

The anti-tagged jet purity Pg is defined to be the fraction of anti-tagged jets that are 
indeed associated with an underlying gluon jet, and the tagging rate, Vtag, is defined to 
be the fraction of normal-mixture Y-events which contain jets that are anti-tagged. From 
Monte Carlo studies, the normal-mixture sample of jets had a gluon content, pn.mix = 
48.7 ± 0.2%. The anti-tagged jet purity and the tag rate depend on the values of 6 and 
/q"*, and values of pg of up to about 80% can be achieved whilst maintaining a tag rate in 
excess of about 30%. A high gluon purity is desirable for the algebraic correction procedure 
later, and therefore G = 7° and /g,"* = 0.75 were chosen for the standard tag. From 
studies of Monte Carlo events the values Ptag = 30.3 ± 0.2% and Pg = 77.0 ± 0.5 ± 0.9% 
were determined, where the first error is statistical and the second systematic. The 
systematic error includes contributions from the choice of Monte Carlo model (Jetset 7.4, 
Jetset 7.3 or Herwig 5.9) and the method used to identify the quark jet in the simulated 
events (the method described in section p.2| for the energy-based analysis was used as an 
alternative). In the data, 6 = 7° and f^""' = 0.75 gave a tag rate of 31.0 ± 0.2% (23,256 
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anti-tagged jets) which is well described by the simulated events. The average energy of 
the anti-tagged jets was 23.59 ± 0.02 GeV. 

In about 1.9% of events both of the lower energy jets fulfilled the tag criteria due to 
misidentification of collimated gluon jets (the highest energy jet is the gluon in fewer than 
10% of these cases), and therefore in these events both the lower energy jets were included 
in the anti-tagged gluon jet sample. The double tag rate rose to about 3.5% for 6 = 8° 
and /I,"* =0.75. 

The flavour composition of the tagged, anti-tagged and normal-mixture jets in the 
Monte Carlo events is shown in table ^. Significantly fewer b-flavour jets were tagged 
than light-flavour jets which is consistent with the observation that b quark jets in Z° 
decays are broader than light quark jets [Q; a small reduction is also visible in the rate of 
tagging c-fiavour jets. The light quark flavour composition of the tagged jets reflects the 
relative couplings of u- and d-type quarks to the Z°. There is a fairly good correspondence 
between the quark flavour properties of the normal-mixture and anti-tagged jet samples. 
It was shown in that there is no significant systematic bias to the measurements of 
quark and gluon jet differences from the flavour composition of the anti-tagged jets. The 
flavour composition of the quark jets in the anti-tagged sample is therefore not expected 
to give a systematic bias in the present analysis. 



4.3 Correction Methods 
Strange Particle Identification 



Kg mesons and A baryons were reconstructed as described in section |2]^. Their invariant 
mass spectra were computed in bins of momentum, andpA, within the normal-mixture 
and anti-tagged jets separately. A Kg or a A candidate was assigned to a jet if its 
momentum vector fell within the cone defining the jet, i.e., if it lay less than 30° from the 
jet axis. If a Kg or A was within 30° of both jets, it was assigned to the higher (calculated) 
energy jet. 

The detection efficiencies of the Kg and A finding algorithms were sensitive to the 
particle momenta and to the track environment, and were therefore determined from 
simulated events for each bin of p^o and p\ within the normal-mixture and anti-tagged 
jets separately. The efficiency was computed from the fraction of generated Kg or A 
within a jet that were found in the same jet by the algorithm. The background-subtracted 
numbers of Kg and A were corrected for detection efficiency as a function of momentum, 
and summed to give the production rates of Kg and A in normal-mixture and anti-tagged 

jets, -Dn.mix ^'^d -Da.tag; -Dn.mix ^'^cl -Da.tag- 



Algebraic Correction Procedure 

The algebraic method introduced in P| was used to correct for quark and gluon misiden- 
tification and to arrive at the ratio of the identified particle production rates for pure 
quark and gluon jets. 
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The production rate of particle type i in the normal-mixture sample of jets, -Dn mix 
may be written 

-^n.mix = Pn.mix ' + (1 — Pn.mix) ' , (2) 

where and are the production rates of i in pure gluon and quark jets respectively. 
Similarly, the production rate in the anti-tagged sample -D^ ^^^g may be written 

^Itag = Pg ■ + (1 - Pg) ■ Q' . (3) 

The ratio of the production rates for pure gluon and quark jets 71^^ may therefore be 
determined by 

_ G' _ (1 - Pn.mix) ■ (I^ltag/^n.mix) " (1 " Pg) 
^"^ Q' Pg - Pn.mix (-DLtag/-D 

n.mix/ 

As in the case of the multiplicity measurement in |^ it was assumed for the purposes 
of the algebraic correction that and are the same in equations |^ and ^ This is 
a reasonable assumption, as the quark flavour compositions of the normal-mixture and 
anti-tagged samples are generally consistent, and as the properties of energetic acolinear 
gluons are independent of event flavour according to QCD. The simulated events provide 
a good representation of the relevant event properties such as the coUimation of jets. 
Therefore any residual effects are expected to be removed by the corrections for detector 
effects. 

The measured production rates appear in equation ^ as the ratio (^'ItagZ-^mmix) 
therefore some systematic effects related to particle identification are expected to cancel. 
Statistical uncertainties on TZ^^^ were estimated from the variance of the results obtained 
when the analysis was repeated ten times with the data split into independent subsets. 
This procedure correctly takes into account correlations between Dl ^^^ and -Dn. mix- 
Detector Corrections 

A correction derived from Monte Carlo events was applied to correct for detector accep- 
tance and resolution. The correction was formed from the ratio of TZ^^ values determined 
at the generator and detector levels. At the generator level, the same three-jet event 
selection criteria as for the data were applied (with the exception of the requirement 
|cos^| < 0.9 for the jet axes). Monte Carlo information was used to identify quark 
and gluon jets, as well as particle type. The detector corrections were determined to be 
1.105 ± 0.005 for charged particles, 1.023 ± 0.006 for Kg and 1.058 ± 0.008 for A with 
uncertainties due to the limited statistics of the Monte Carlo samples. The correction 
factor for charged particles agrees well with that determined in 0. The final results can 
be found in table ^ which will be discussed later. 



4.4 Systematic Uncertainties 

Systematic uncertainties, which are listed in table H, were evaluated from a number of 
sources in a similar fashion to the energy-based analysis described in section The total 
systematic errors were obtained by combining the individual contributions in quadrature. 
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Detector effects: The total energy and momentum flow in the event were estimated 
using charged tracks only instead of tracks and electromagnetic clusters. These 
were used as input to the jet finding and quark jet identification algorithms and 
the analysis was repeated. The observed jets are more coUimated when measured 
using tracking information only. This is reproduced by the Monte Carlo, both at 
the generator and detector levels and is partly due to decays of 7r° into photons 
which tend to broaden the flow of neutral particles in the jet. The effect of this 
increased collimation was that the anti-tagged jet purities and tag rates were some- 
what different to the standard case, with pg = 70.2 ± 0.4% and Ptag = 42.0 ± 0.2% 
for = 7° and /q"* = 0.75. The uncertainty was estimated from the difference 
divided by vT2/2 as for the energy-based analysis above. 

Quark jet identification: The analysis was repeated using parameters G = 8° and 
/e"* = 0.75 since its tagging rate ("Ptag = 40.1 ± 0.2), and the anti-tagged jet purity 
(pg = 74.7 ±0.4 ±0.9) were somewhat different compared to the standard analysis. 

Event selection: Events were selected requiring a minimum jet energy of 10 GeV which 
is the selection criterion used in 

Jet purity determination: A number of sources of systematic uncertainty on pg have 
been considered. Simulated events from the Jetset 7.4 and Jetset 7.3 samples 
were used separately to determine the purities and tag rates, and the observed dif- 
ferences were used to evaluate an uncertainty related to the event generator tuning 
and detector simulation. A further uncertainty comes from the ambiguity in defin- 
ing whether a jet arises from a quark or a gluon, as described in a previous OPAL 
publication Q . The same procedures were followed to determine a systematic error 
which was added in quadrature with the other sources to arrive at the total system- 
atic uncertainty. The analysis was repeated using the tagged and normal-mixture 
jet purity values varied by their combined statistical and systematic uncertainties. 

Background determination: The fit ranges, signal window sizes and excluded regions 
were all varied to determine the background to the selected Kg and A signals. In 
the case of the A an alternative function was also used to describe the background. 
Finally a sideband method was used to determine the background. 

Efficiency determination: The Jetset 7.3 and Jetset 7.4 samples of simulated events 
were used separately to estimate the efficiency of the Kg and A finding algorithms. 
The two Monte Carlo samples are taken to represent alternative possibilities of 
generator tuning and detector simulation and an uncertainty was estimated as for 
the energy-based analysis above. Several of the systematic variations were made 
simultaneously, and in no case was a difference larger than that for Jetset 7.3 
observed. 

Monte Carlo statistics: The finite numbers of simulated events available led to statis- 
tical uncertainties in the particle detection efficiencies, and the detector corrections. 
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The effect of varying the cone size R that defines the jets has been investigated in |^ 
where an increase in T^g^ with R was observed. Other sources of uncertainty such as 
changing the requirement on the angle between the highest energy and the other jets, and 
modifying the track and cluster selection criteria have been considered in that publication 
and found to be negligible. There were also no statistically significant differences between 
the detector acceptance corrections computed with the Jetset 7.3 or Jetset 7.4 samples. 



5 Results and Discussion 

The results of the energy-based analysis are given in table The ratios Rg^/Rq^ and 
Rg/Rg of relative Kg and A production rates in pure gluon and quark jets are given for 
all three event selections, together with their statistical and systematic uncertainties and 
the predictions of the Jetset 7.4 model. For the |/cut selection, the production of Kg 
mesons is enhanced in gluon jets relative to quark jets by a factor 1.10 ±0.02 ±0.02 where 
the first error is statistical and the second systematic. Similar results were obtained with 
the other event selections. The relative production rate of A baryons was found to be 
increased in gluon jets relative to quark jets by 1.41 ± 0.04 ± 0.04 in the ycut selection. 
The less collimated y„in and cone event samples show a smaller increase in the relative 
production rates of A baryons indicating a possible dependence of the baryon production 
rates on the topology of the events. 

For the Y-event analysis, the ratios of absolute production rates of Kg and A in 24 GeV 
gluon and quark jets were found to be 7^g = 1.05±0.08±0.09 and 7^^q = 1.32±0.11±0.18 
for jets defined with the cone algorithm using a cone size of 30°. The corresponding result 
for charged particles is T^gq = 1.116 ± 0.006 ± 0.012 which is in good agreement with the 
result from the previous OPAL vertex-tagged analysis (1.10 ± 0.02 ± 0.02 [Q) and which 
has a reduced statistical error as a result of the use of the more efficient energy flow tag. 

The measurements in the Y-event analysis may also be used to determine production 
rates of Kg and A in quark and gluon jets relative to those of charged particles. The 

relative rates of K° and A production are i?g Vi?P = 0.94 ± 0.07 ± 0.07 and R^/Rf = 
1.18±0.10±0.17, where correlations between the sources of systematic uncertainty (given 
in table ^ have been taken into account. 

An enhancement of A production in gluon jets relative to quark jets, in excess of that 
observed for charged particles, is measured by both analyses, with the ratios of the relative 
A rates consistent within the errors. The ratios of the relative production rates of Kg 
mesons in gluon and quark jets are also compatible within the errors, and suggest a small 
enhancement relative to charged particles. The analyses presented here are consistent 
with previous findings that have shown that factors such as jet finder, jet energy 



and event topology are important in quantifying the differences between quark and gluon 
jets. Therefore, care should be taken to ensure that the conditions are equivalent when 
comparing results between experiments. 

The measurements of relative production rates in gluon and quark jets from both 
analyses are shown in figure |^ together with the predictions of the Jetset 7.4 and Her- 



18 



WIG 5.9 Monte Carlo models. These data are also given in table ^ together with the 
measurements of absolute rates from the Y-event analysis. Herwig 5.9, despite its good 
description of global event properties, fails to give an adequate description of the mea- 
surements of ratios of strange particle production rates in quark and gluon jets. This 
result is not surprising given the poor description of the inclusive strange particle rates as 
a function of jet energy shown in figure |^. Jetset 7.4, however, was shown to provide a 
reasonable general description of the data in figure and the ratios of relative rates from 
the Y-event sample are consistent with its predictions for both Kg and A. The ratios of 
relative rates from the energy-based analysis, however, are significantly larger than those 
predicted for both Kg and A. The predictions of the ratios of relative production rates 
from Jetset 7.3 differ from Jetset 7.4 by at most 0.04 which gives an indication of the 
size of possible effects due to parameter tuning and inclusion of additional particle decay 
channels. 

There is no perturbative mechanism in the Jetset model that gives rise to the ob- 
served differences in particle production between quark and gluon jets; they arise from 
the effects of hadronization and particle decays. Many of the Kg are the decay products 
of heavy (b and c flavour) hadrons. As the production of b and c quarks in gluon jets 
is suppressed, there are correspondingly fewer Kg in these jets. The enhancement of A 
production (and of baryons in general) in gluon jets relative to quark jets predicted by 
Jetset is a consequence of the different dynamics of the string fragmentation process in 
quark and gluon jets. 

6 Summary 

Production rates for Kg and A have been measured in quark and gluon jets from Z° decays 
with two complementary approaches. In the first analysis a procedure was introduced to 
compare particle production in gluon and quark jets of different energies facilitating the 
study of up to about 24% of the total event sample. Different jet finding algorithms {k±, 
cone) were used to investigate samples with different three-jet event topologies. Relative 
rates, normalized to the inclusive charged particle rate, were obtained for pure quark 
and gluon jets considering the different quark and gluon content of jets in different energy 
intervals. The relative rates in pure gluon and quark jets were found to be consistent with 
being independent of the jet energies, and to depend slightly on the specific jet selection. 

In the second analysis a new method was introduced to tag quark jets based on the 
collimation of their energy flow, allowing the isolation of a larger sample of anti-tagged 
gluon jets in symmetric three-jet events than the method of secondary vertex tagging used 
previously. The comparison of quark and gluon jets of equal energies and embedded in 
almost identical event environments allows for a simple interpretation of the results. The 
jets were selected using a cone algorithm. By also measuring the inclusive particle rates 
in these symmetric jets, relative rates were obtained in addition to the absolute rates. 

An enhancement of A production in gluon jets relative to quark jets, in excess of that 
for charged particles is observed. The ratios of production rates of Kg mesons in gluon and 
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quark jets suggest a small enhancement relative to charged particles. The results of both 
analyses are compatible within their errors. The predictions of Jetset 7.4 are consistent 
with the enhancement observed for A, but are smaller for the Kg, whilst Herwig fails to 
provide an adequate description of the data. 
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Tables 



Energy-based 


Kg 


rLq 


xVd.o.f. 


Vcut 


0.0573 ±0.0009 


0.0522 ±0.0006 


9/13 


Z/win 


0.0568 ±0.0009 


0.0531 ± 0.0006 


8/13 


cone 


0.0622 ±0.0011 


0.0580 ±0.0007 


14/13 


Energy-based 






xVd.o.f. 


Vcut 


0.0252 ±0.0005 


0.0179 ±0.0003 


15/13 


I/win 


0.0252 ±0.0005 


0.0188 ±0.0003 


12/13 


cone 


0.0281 ± 0.0006 


0.0214 ±0.0004 


39/13 



Table 1: Fitted relative Kg and A production rates in gluon and quark jets from the 
different event selections of the energy-based analysis. The results have been obtained 
by a fit to equation 1; the fit quality is indicated by the x^/d.o.f. values. The errors are 
statistical. 



Energy-based 


Rg^/Rq^ 


Jetset 7.4 


Z/cut 


1.10 ±0.02 ±0.02 


0.94 




1.07 ±0.02 ±0.02 


0.94 


cone 


1.07 ±0.02 ±0.02 


0.95 


Energy-based 




Jetset 7.4 


ycut 


1.41 ±0.04 ±0.04 


1.26 


Z/win 


1.34 ±0.04 ±0.03 


1.24 


cone 


1.31 ±0.04 ±0.05 


1.30 



Table 2: Ratios of the relative Kg and A production rates in gluon and quark jets from 
the different event selections of the energy-based analysis compared to the Jetset 7.4 
predictions. The first error is statistical and the second systematic. 
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Source of error 




A 


Statistical error 


2.1 % 2.3 % 2.3 % 


2.7 % 3.0 % 3.2 % 


Charged tracks only 
1/cut/l/win change 


0.1 % 
0.2 % 


0.4 % 
0.2 % 


A/Kg cut variation 
Sideband fit 
JETSET 7.3/7.4 


1.1 % 

0.6 % 
0.4 % 


1.7% 

0.1 % 
0.9 % 


Fit range 
Jet purities 
Fit method 
MC slopes 


0.7 % 0.7 % 0.9 % 
0.3 % 0.3 % 0.4 % 
0.6% 0.7% 0.5% 
0.3% 0.2% 0.8% 


1.4% 0.5% 1.4% 
0.8 % 0.4 % 0.5 % 
0.7 % 0.8 % 2.2 % 
0.7% 0.7% 1.2% 


Total systematic error 


1.7% 1.8% 2.0% 


2.7 % 2.3 % 3.5 % 



Table 3: Statistical and systematic errors of the ratios of relative particle production rates 
in gluon and quark jets in the energy-based analysis. The three numbers in a row refer 
to the ^cut) Uvrin, ^ud cone selection. 





Anti-tagged 


Tagged 


Normal- mixture 


d 


5.1 ±0.1 


22.1 ±0.2 


11.3 ±0.1 


u 


4.2 ±0.1 


18.4 ±0.2 


9.1 ±0.1 


s 


5.2 ±0.1 


23.4 ±0.2 


11.4 ±0.1 


c 


4.2 ±0.1 


12.8 ±0.1 


9.0 ±0.1 


b 


4.2 ±0.1 


4.3 ±0.1 


10.5 ±0.1 


gluon 


77.0 ± 0.5 


19.0 ±0.2 


48.7 ±0.2 



Table 4: Flavour composition (in %), in the Y-event analysis, of the anti-tagged, tagged 
and normal- mixture jet samples determined from the Jetset Monte Carlo including 
simulation of the detector. The quark jet tagging used the parameters 6 = 7° and /q"* = 
0.75. The errors are statistical only. 
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Y-events 


Absolute rates 


Relative rates 




-pch 








A 


Detector effects 


0.003 


< 0.01 


0.04 


< 0.01 


0.04 


Quark jet identification 


0.002 


0.05 


0.07 


0.04 


0.06 


Event selection 


0.010 


0.03 


0.08 


0.02 


0.08 


Jet purity determination 


0.003 


< 0.01 


0.03 


< 0.01 


0.01 


Background determination 




0.04 


0.05 


0.04 


0.05 


Efficiency determination 


< 0.001 


0.04 


0.11 


0.04 


0.10 


Monte Carlo statistics 


0.005 


0.03 


0.07 


0.03 


0.07 


Total systematic error 


0.012 


0.09 


0.18 


0.07 


0.17 


Statistical error 


0.006 


0.08 


0.11 


0.07 


0.10 



Tabic 5: Breakdown of the contributions to the uncertainties on the ratios of particle 
production in quark and gluon jets from the Y-event analysis. 



Ratios of relative rates 


OPAL Data 


Jetset 7.4 


Herwig 5.9 


Energy-based 




1.10 ±0.02 ±0.02 


0.94 


0.73 


iVcut) 


A 


1.41 ±0.04 ±0.04 


1.26 


0.88 


Y-events 




0.94 ±0.07 ±0.07 


0.95 


0.62 




A 


1.18 ±0.10 ±0.17 


1.34 


0.87 


Ratios of absolute rates 


OPAL Data 


Jetset 7.4 


Herwig 5.9 


Y-events 




1.05 ±0.08 ±0.08 


1.06 


0.70 




A 


1.32 ±0.11 ±0.18 


1.56 


0.99 




charged 


1.116 ±0.006 ±0.012 


1.16 


1.13 



Table 6: Ratios of the relative Kg and A production rates in gluon and quark jets from 
both analyses, together with the predictions of Jetset 7.4 (the predictions of Jetset 7.3 
correspond to those of Jetset 7.4 for Kg and are about 0.04 lower for A) and Herwig 5.9. 
The ratios of absolute rates determined in the Y-events analysis are also shown. 
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Figure 1: Jet energy distributions of thrcc-jct events selected with the k± jet finder with 
the ^cut event selection. The jets are ordered according to their assigned energies. 
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jet energy [GeV] 



Figure 2: Purities of the reconstructed jets in Monte Carlo events for the i/cut event sample. 
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Figure 3: Production rates per jet of charged particles, Kg mesons, and A baryons, Wc/i/lO, 
rij^Q, and ua, from the t/cut sample as a function of the jet energy compared with the 
predictions of the models Jetset 7.4 and Herwig 5.9. The charged particle rates are 
scaled down by a factor of 10. The errors shown are the (uncorrelated) statistical ones 
and are mostly smaller than the size of the symbols. 
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Figure 4: Relative production rates of (a) Kg and (b) A in Jetset 7.4 events for pure 
quark and gluon jets as a function of the jet energy. The lines are fits of straight lines to 
the points. The statistical errors are smaller than the size of the symbols. The zeros of 
the vertical axes have been suppressed. 



29 



o in 

^ s 



0.065 



0.06 



0.055 



0.05 



0.045 



-i — I — I — I I I I — |- 

a) Ks 



T — I — I — I — I — r 



I ' ' ' ' I 



T — I — I — I — r 



I ' ' ' ' I 



OPAL 




■ OPAL data 
— fitted function 



J I'll 1 1 1 1—1 I I I I I I'll 1 1 1 1 1 1 I I I I 1 1 1 I'll 



10 15 20 25 30 35 40 45 50 

jet energy [GeV] 



< s 



0.03 



0.025 



0.02 



0.015 



0.01 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

: b)A i 

1 

' T— — • 


1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

OPAL \ 


• OPAL data 


— 

i — 1 
• 


1 

L • 


— fitted function 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



5 10 15 20 25 30 35 40 45 50 

jet energy [GeV] 



Figure 5: Relative production rates of (a) Kg and (b) A from the ?/cut selection as a function 
of the jet energy. The lines show the functions returned from the fits of equation |I] to the 
data. The zeros of the vertical axes have been suppressed, and the errors shown include 
both statistical and systematic uncertainties. 
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Figure 6: Distributions of /e for quark and gluon jets. The events were generated with 
the Jetset Monte Carlo and include full simulation of the detector. The arrow indicates 
the cut position for selecting quark jets. 
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Figure 7: The ratio of relative production rates (see text) in quark and gluon jets of 
Kg and A for both analyses. The experimental statistical errors are delimited by the 
small vertical bars. The predictions of Jetset 7.4 and Herwig 5.9 are also shown. The 
predictions of Jetset 7.3 are no more than 0.04 lower than those of Jetset 7.4. 
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